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We show that two commonly used definitions for the heat current give different resuhs — through 
the Kubo formula — for the heat conductivity of oscillator chains. The difference exists for finite 
' chains, and is expected to be important more generally for small structures. For a chain of A'' 

, particles that are tethered at the ends, the ratio of the heat conductivities calculated with the two 

• currents differs from unity by 0{1/N). For a chain held at constant pressure, the difference from 

^S) ' unity decays more slowly, and is consistent with 0{1/N^) with 1 > > 0.5. 
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I. INTRODUCTION 



In linear response theory, transport currents that flow in a system in response to small gradients in thermodynamic 
potentials are calculated in terms of equilibrium autocorrelation functions of these currents through the Green-Kubo 
^ formula This involves taking the thermodynamic limit of the autocorrelation functions first, and then the zero 
D frequency limit; the order in which the limits are taken is important The response function to be calculated 
^ I determines the appropriate correlation function: for instance, the electric and thermal conductivities involve the 
^ I ^ ' autocorrelation functions of the particle and heat currents respectively. 

. Unlike the particle current, which is defined unambiguously, there are many possible choices for the energy current, 
' and therefore the heat current which is a linear combination of the two Q. Most notably, the heat current can be 
^ ; , defined in the continuum or on a lattice; both are commonly used. For the lattice current, the energy of each particle 
d is assigned to the lattice site associated with the particle. This is commonly used for oscillator chains, and seems 
^ \ reasonable for crystals (if we neglect lattice defects). The energy current fiowing from the lattice site i to A: is 

■ j»fc = iF,fc-(v, +Vfc) (1) 

c ■ 

O , where Fik is the force exerted on the fc'th particle by the i'th particle, and v^^fc are the velocities of the particles. For 
the continuum current, commonly used for hard particle systems or fluids, the energy of each particle resides at its 
instantaneous location. As a result, there is an advective part to the energy current from the motion of particles 



^! jadt,(x) = ^e.iVi(5(x- Xi) (2) 

o 

' where is the energy of the i'th particle. In addition, since the energy current, jit, between particles i and k flows 

. 1 between Xk and Xi instead of between their lattice sites, the spatial integral of jik (used in the Green-Kubo formula) 

' is also different. 



In view of the obvious approximations in the lattice current, is it merely something that works well for stiff 
crystals? When computing the conductivity from equilibrium correlation functions, are the two currents equivalent in 
the thermodynamic limit [5|? This would not be sufficient for small structures that are s pec ially designed and cannot 
be scaled up, for which the conductance has to be used instead of the conductivity [l6l. Il7||. Even for large systems, 
k>i ^ one has to be careful in proving the equivalence of the two currents if the transport coefficients are singular, since the 
5_j thermodynamic limit of the conductivity does not exist. 

_ Surprisingly, despite the approximations in the lattice current, one can prove an exact Green-Kubo like formula 
with it for the conductivity of a finite chain which has Langevin baths attached to the terminal particles Q , using the 
steady state fluctuation theorem [l^. This result has been proved without the fluctuation theorem and generalized 
to a number of different implementations of heat baths [Toj . The proof applies to a chain with arbitrary onsite and 
interparticle potentials (which may vary down the chain, even though the proof does not state this). Therefore, it 
applies to a tethered or free chain and (by a simple transformation that we shall show in this paper) a chain at 
constant pressure. Is this agreement accidental, and can the proof of Ref. 0] be extended to the (exact) continuum 
current as well? 

In this paper, we explain the result by obtaining an interpretation of the dynamics of the oscillator chain for which 
the lattice current is exact; in this interpretation, the Langevin baths are at fixed locations. We further show through 
analytical and numerical calculations on such a chain that the conductivities obtained using the lattice and continuum 
currents are not equivalent except in the thermodynamic limit. The ratio of the two differs from unity by ~ 1/N for 
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a chain with N particles that is tethered at both ends. For a chain maintained at constant pressure, the approach to 
unity is slower. If fit to a ^ form, it is consistent with any i < 77 < 1. Our numerics also show that slight changes 

in the definition of the continuum heat current, which is a linear combination of the energy and particle currents, can 
markedly change the error in the resultant conductivity; these differences are inconsequential in the thermodynamic 
limit. The conductivities obtained from lattice and continuum currents are generally very close for N > 256. 

One can try to explain the error in the Green-Kubo formula using the continuum current by noting that the size 
of the system fluctuates in the continuum interpretation. In the tethered case, the terminal particles are adjacent to 
the tethering posts and the percentage fluctuations in the length are 1/-/V, while for a chain at constant pressure 
they are ~ 1/\/N. This expectation is only partly borne out by the numerics. Note that, in view of the proof in 
Refs. it is the lattice current that gives the heat conductivity accurately. 

The rest of this paper is organized as follows. Section II reviews the derivation of the lattice and continuum heat 
currents and shows how the lattice heat current is not an approximation if the oscillator dynamics are interpreted 
suitably. Section III calculates the ratio of the conductivities obtained with the continuum and lattice currents and 
how they behave as the thermodynamic limit is approached. Section IV presents numerical results with the same for 
a tethered chain and a chain at constant pressure. 

II. TRANSVERSE AND LONGITUDINAL INTERPRETATIONS 

We consider a one dimensional chain of particles with nonlinear springs connecting adjacent particles. (The defini- 
tions in this section are standard but we include them for clarity.) The Hamiltonian is 

N N-1 N 

i—l i—1 i—1 

where mi,Xi and Vi are the mass, position and velocity of the i'th particle and U is the interparticle potential. Uq is 
an onsite potential for each particle. When Uq = 0, this is a Fermi-Pasta-Ulam (FPU) chain 12]. The dynamics of 
this system are Hamiltonian, with heat baths connected to the first and last particle through the addition of Langevin 
damping — 7Wi,Ar and noise 771, at (t) terms in their equations of motion. As with other one-dimensional systems in 
which interparticle interactions are momentum conserving [lo| , the FPU chain has a singular heat conductivity, 
with the thermal conductivity of the chain growing as ~ N^^^ as a function of chain length jl3, [H, [3j ■ In higher 
dimensions or if Uq 7^ 0, the singularity in the heat conductivity is eliminated. 

If we associate the kinetic energy of each particle with itself, and the potential energy of each spring as being equally 
distributed between its two neighbors, the rate of change of the energy of particle i 

Ci = ^mvi + Uoixi) + ^[U{xi+i - Xi) + U{xi-i - x^)] (4) 

is 

ei = niiViVi + UQ{xi)vi + ^Fis-i{vi-i - Vi) + ^Fis+i{v,+i - v^) (5) 

where i^i.i-i is the force on the z'th particle from the i ~ I'th particle. (There is only one F term in ei and cat.) 
Using the equations of motion, this is equal to 

ei = ^Fis^i{vi^i + Vt) + lFi^i+i{vi+i + Vt.) + Si^i{-^vi + t]i)vi + 5i^Ni-lVN + r]N)vN (6) 

Therefore the continuum energy current is 

N N 

j|(x,t) = ^e^ViS{x - Xi{t)) + ^ ^F,;,i-i(wi-i + Vi)e{x - Xi^i)e{xi - x) (7) 

i=l i=2 

where the first term comes from the advective transport of the energy associated with each particle by it. Integrating, 
the energy current flowing through the entire system is 

N N 

Jsii) = '^^iVi + \ '^Fi,i^i{vi-i +Vi){xi -Xi-i). (8) 

i=l i=2 
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With other choices of how the potential energy of each spring is distributed, this expression can change by small 
amounts. This result can easily be generalized to higher dimensions: 

N 

JfiW = ^e,v, + i^^F,fe •(vfe+v,)(x, -Xfc). (9) 

On the other hand, if the springs between the particles are very stiff, each particle only deviates slightly from its 
lattice position. If we neglect the advective term in the heat current and approximate Xi « ia (where a is the lattice 
constant), we obtain an energy current equal to Fi^i-i(vi-i + Vi) between the i ~ I'th and i'th lattice sites, which 
integrates to the lattice energy current 

N 

4W = ^^^^M-i(^^«-i +^^0a. (10) 

i=1 

The generalization to higher dimensions is again straightforward. For a one-dimensional system, in which the ordering 
of particles is fixed, it is always possible to view the positions of the particles as displacements from a reference lattice, 
but in higher dimensions this is only useful in a crystalline phase. 

By its construction, the lattice energy current seems clearly an approximation; if a Green-Kubo like formula 
exists for finite oscillator chains, one would expect it to be only approximate if the lattice current is used. Despite 
this, one can prove an exact Green-Kubo type formula for finite oscillator chains with Langevin baths if the lattice 
current is used 

^=^^^;,(^/°°(^g(t)JQ(0))d. (11) 

In Eq. pT|) Jq(€) is the spatially integrated heat current (defined in the next section) flowing at time t through the 
chain of N particles at a temperature T. 

As mentioned in the previous section, one can construct an interpretation of the oscillator chain dynamics for which 
the lattice current is exact. Let Xi = ia + where Q is the displacement of the i'th particle from its lattice position. 
We imagine that the displacement is transverse to the chain instead of along it. If the potential between neighboring 
particles is expressed as U{Ci+i -fa — Q), the dynamics are the same as if the displacements Ci had been in the 
longitudinal direction. The energy current flowing into the Langevin baths is also unchanged. In this interpretation, 
the spatially integrated energy current is with no approximations (20| : there is no advective term in the current, 
and the energy current between neighboring particles flows over a distance a instead of Xi — Xi-i = a + Q — Ci-i • (The 
transverse interpretation is in fact more natural for the often-studied 'phantom' oscillator chains, where particles can 
pass through each other but only neighboring lattice sites interact.) 

Since the lattice and continuum currents are exact for the transverse and longitudinal interpretations respectively, 
and the heat flowing between the reservoirs is independent of interpretation, one might hope that Eg. ljlip would be 
valid for both currents. In the rest of this paper, we shall show that this is not the case except in the thermodynamic 
limit. 



III. LATTICE AND CONTINUUM HEAT CURRENTS 

The proof in Ref. uses the lattice energy current for the heat current for a system at zero pressure. Applying 
a pressure at the ends of the chain is equivalent to replacing U{xi — Xi^i) with U{xi — Xi-i) — p{xi — Xi^i). From 
Eg. pO)) . this changes the heat current to 



N 

J^t) = lEt^M-i (12) 

i=2 

with which definition the proof of Ref. [6l| is extended to systems at non-zero pressure. The heat current in the 
continuum interpretation is defined through Galilean invariance as 0, 0] 

Jqit) = Jsii) -[E+{N~ l)pa\vcM (13) 

where E is the time average of the energy of the system and vcm is the instantaneous velocity of the center of mass 
defined by 



1 ^ 1 



N N 

i=l 
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where Jjv is the number current. 

In comparing Jg and Jq, it is useful to first derive an identity that is vahd for the autocorrelation function of any 
conserved current: 



1 

C{t)dt = lim — 



^ '■^ 1 / r 1 

J{ti)-J{t2))dtidt2^ lim -( I 3{t)dt 



Jo 







(15) 



Thus in Eq.pT]). we are interested in the 0(\/t) P^rt of 3Q{t)dt. In particular, in Eg. p^ we have 

N 

lP^i'"i~i + Vi)a = {N - l)pavcM 



(16) 



i=2 



similar to Eg. p^ . since /J"[i'i(i) — vcM{t)]dt is equal to the change in Xi{t) — xcM{t), which cannot be 0{^/t). Note 
that even though Jqit) has a power law tail to its autocorrelation function, this tail is cut off for large time for any 
finite TV, as it must be from Eq. p^ if the thermal conductivity is to be finite for that N . 
As another useful result, if p is the charge density corresponding to a current jp, then 



d 
di 



■ J xp(x, t)dx = - y xV • jp(x, i)dx = - J xjp(x, t)-dS + J pit) (17) 
in which the integral on the right hand side is over the surface of the system, where it is connected to reservoirs. 



A. Tethered chain 



Let the FPU chain of Eq.® be tethered at the ends. This is accomplished by adding extra fixed particles at the 
zeroth and TV + I'th locations with springs connecting them to their neighbors. With these boundary conditions, the 
displacement of the center of mass in a time interval r — + oo must be 0(1). Then from Eg.psp. the energy and heat 
currents yield the same result in Eg. pip , both for the lattice and continuum currents. 

From the continuity equation, 

d ^ AN-l)a ,-i , ,N 

-Y^{^-l)ae,{t) = - (y ^ a)^^^dy ^ f^it) - {N - l)aJEAt) (18) 

where j/ is a continuous coordinate along the lattice, such that j^(y) is piecewise constant with discontinuities at 
integer multiples of the lattice constant a. Js-R it) is the energy current flowing into the heat bath at the right end of 
the system. If we integrate both sides of this equation over an extremely large time interval r, the left hand side does 
not diverge with r. The first term on the right hand side must be 0{y/T) from Eq. lTS]) . and therefore so must be the 
second term. Keeping only terms that grow as 0{^/t), we obtain 

r />r 

fE{t)dt ={N-- l)a / JE;R{t)di- (19) 

Similarly, since from Eq.lISl) the left hand side of Eq.^H) is 0(TV(i+")/2) for large N where a is the heat conductivity 
exponent [ll|, the right hand side must be the same. 

We can try to understand the 0(-y/rA^^^+"^/^) scaling of the right hand side as follows. On extremely long time 
scales, J_E;ij(t) is as likely to be negative as positive. Thus the right hand side of Eg. ([TO)) is 0{^/t) for large r. On 
the other hand, if energy flows into the system from the reservoir to the right, it increases the local energy density, 
thereby increasing the likelihood that energy will flow out shortly afterwords. As TV ^ oo, all the energy that flows 
in from the reservoir to the right must eventually flow out to the same reservoir instead of escaping to the left. 
Consequently, Js-Rit) is anticorrelated up to a time scale that diverges with TV. This reduces the 0(TV) dependence 
one might naively expect for the right hand side of Eq. (fTQ]) . Although the arguments in this paragraph are qualitative, 
the 0(-yTTV(^+")/^) scaling of the right hand side of Eg.p^ is derived rigorously in the previous paragraph. 

For the continuum energy current, 

d x-^ d f^^^^'i f dj'^ (x t) 

— > ei(t)a;^ = — / e{x,t)xdx = - x^^^j-^dx ^ Jl^it) + xi{t)JE;L{t) - XNit)JE:Rit) (20) 
dt ^ dt J^^i^t^ J dx 
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where J el is the heat current flowing out of the reservoir on the left hand side of the system. With xi N{t) = 
[0,(iV-l)a]+Ci,A'(i) 

' rE{t)dt ={N^ l)a r JEMt)dt + I [CiV " Cl(t)JE;L{t)]dt. (21) 



Compared to Eg. p^ . the extra term on the right hand side is also 0{\/t) for large r. However, the missing factor of 
[N — 1) makes it negligible compared to the first term in the thermodynamic limit. 

One would expect the factors of C, in the integrand to destroy the temporal anticorrelation of Je;L and Je-.r so that 
the last term should be independent of N. Thus we would expect 

/o^^QqW^^ _ ^ _ _ ^-4/3^ (22) 

Contrary to this expectation, we will see numerically that the right hand side is in fact proportional to as 
suggested in Section|ll We must conclude that the correlations between ^ and Je;L,r are more subtle than one would 
naively expect. 



B. Constant pressure boundary conditions 



If the FPU chain has a constant pressure p applied to it instead of being tethered at the ends, the equivalence of 
the continuum and lattice heat currents is more delicate. The tethering is removed, and the pressure effectively adds 
a term p{xpf — xi) to the potential energy in Eq.Q. The center of mass of the system executes a random walk due 
to the fluctuating forces exerted by the heat baths, moving a distance 0(-\/r) in time r. The volume of the system 
also fluctuates by 0{\/N). However, the lattice description is essentially the same as for the tethered chain. Since 
the system at a pressure p is equivalent to a system at zero pressure with an extra p{xi — Xi-i) in the interparticle 
energy, from Ea. (jl2p the continuity equation Eq. (|18p is modified to 

d ^ 

— ^{i - l)a[e^{t) + \pa{x^ - Xr-i) + \pa(x^+l - x^)] = jj5(t) -pa^\{v^+ v^-l) - (TV - l)aJE;R{t) (23) 

4=2 

from which 

jl 



fQ{t)dt = {N ~ l)a JE;R{i)dt (24) 
Jo 

as before. Note that we could have used (i — 1 — k)a on the left hand side of Eg. ([^5]) . which would have yielded a 
linear combination of JE;R.{t) and JE;L{t) on the right hand side of Eq.([M|). 
Turning to the continuum current, 

d / \ Tr/\ T/N ^'^1 dXjV /„ N 

-j^Z^ei{t)Xi = J^it) +xiJE;L(t) -XNJE-RKt) +P^i~j'f ~P^^~dr ^ ' 

where on the right hand side we have included the energy current from the reservoirs and the work done by the applied 
pressure. We define the heat current (integrated over the chain) as 

Jqit) = J%{t) - [E{t) +p[xN{t) - x,{t)]}vcM (26) 

where E{t) is the total energy of the system at time t and vcm is the velocity of the center of mass. This definition 
is slightly different from that in Ea. (fT3)l : if E{t) and XM{t) — xi{t) are replaced with their average values, we recover 
the earlier expression. We also use conservation of energy 

dElt) ^ , , ^ , . dxi dxN ,„„x 

Multiplying this equation by xcm^ subtracting from Eq. (j25p and using Ea. (j26p . we have 
d 



dt 



^ei{t)[xi - Xcm] = Jgit) + {xi - XcM)JE;Lit) - {XN - XcM)JE;R{t) 

+ p{xi- xcm)^{xi- xcm) - p{xN - xcm)^{xn - xcm)- (28) 
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Therefore 

f J^{t)dt = i(7V - l)a ( [JE-Mt) + JE-L{t)]dt + f (Cn- CcM)JEAt)dt + I (CcM - Ci)JE-At)dt. (29) 
Jo Jo Jo Jo 

As per the discussion after Eg. ([24|) . the first term on the right hand side is equivalent to jQ{t)dt. As per the 
discussion after Eg. p^ . this is ~ V tN^+°'. The last two terms are also proportional to y^r. Unlike for the tethered 
case, since Ci,JV — Ccm is 0{^/N), it is not clear they can be neglected even in the thermodynamic limit. Therefore, to 
compare the lattice and continuum conductivities with constant pressure boundary conditions, we turn to numerical 
simulations in the next section. Note that if we had used the first definition of the heat current, Eg. p^ . Eq.(P5|) 
would have had an extra term on the right hand side 

[E - E{t) - pCnit) +pCi{t)]vcM{t)dt (30) 

with which it is even less clear that the lattice and continuum results are equivalent. 



IV. NUMERICAL RESULTS 



In the numerical simulations we set the on-site potential Uo{x) in Eq. ^ to be zero, and take the interparticle 
potential to be 

u{x,+, -x.) = ^ (C.+i - + f (C.+i - CO' + 5 (C.+i - Q)" , (31) 

where Q = Xi — i a The lattice spacing a is taken to be unity. We consider N particles of which the first and last are 
coupled to heat baths at temperature T. Each particle in the interior interacts with neighbors to the left and to the 
right, while the first and last particles only interact with one neighbor. These end particles have additional forces due 
to the friction and noise, — 71^1 + 771 and — 7fAr + ?7Ar, and also interact with an extra tethered 'particle' at the zeroth and 
N + I'th sites respectively. The noise satisfies the fluctuation-dissipation theorem {r]i{t)rjj{t')) — 2kBT^ 5{t — t') Sij. 

We integrated the equations of motion using the second order velocity Verlet (leapfrog) method for all particles in 
the interior. For those on the boundary, which are coupled to the heat baths and so subjected to random noise, we 
just used the simple Euler method, so, for example, for particle 1, 

vi{t + 5t) = v{t) + F(Ci ~ C2)6t -jv{t)6t + Re{t), (32) 

where the force, i^(Ci — C2), is equal to —dU /d(i with U given by Eq. ([3T|) . and e is a Gaussian random variable with 
mean and standard deviation unity. The coefficient of the noise, R, is the root mean square fluctuation in the noise 
integrated over time St, so it is given by 

/•St /'St /'St /'St 

R^ = dt dt' {rj{t)r]{t')) = dt dt' 2kBTj6{t-t') =2kBTj5t. (33) 
Jo Jo Jo Jo 

We focused on one set of parameters, w = 0.5, m = 1.0,7 = 0.5, T = ks = 1.0, and for the constant pressure 
simulations we took p = 0.1. 

We start the system off with all velocities equal to zero, and the particles separated by their equilibrium distance 
a, so all the energy is fed in from the baths at the boundaries. We run for a time icquii to equilibrate, and then 
continue for an additional time t,ncas during which measurements of the current are made after every time step. This 
whole process, both equilibrating and averaging, is then repeated nj-un times and the results averaged. Error bars are 
estimated from the standard deviation of results from different runs in the usual way. The values of the parameters 
used are given in Table |T1 The thermal conductivity is obtained from the simulations from 

where (• • • ) denotes the average over the rij-^n runs. 

We investigated the size-dependence of the difference between the thermal conductivities calculated with the con- 
tinuum currents, Eqs. (fTS]) and (pS)) . from that obtained with the lattice current, Eq. (H^). This difference is always 
small, and, especially for the tethered chain simulations, depends sensitively on the time step, St. We shall therefore 
present results as a function of both St and N. As discussed at the start of Sec. IIII Al for a tethered chain in which 
the heat current is defined by Eq. p3p . the terms subtracted from the energy current are inconsequential for the 
conductance. 
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FIG. 1: (Color online) Results for the tethered chain. A log-log plot of the ratio of the continuum thermal conductivity with 
the heat current defined in Eq. (|13|l to the lattice thermal conductivity with the heat current is defined in Eq. (|12|l minus 
one. If the continuum and lattice thermal conductivities agree in the thermodynamic limit, the results should tend to zero for 
N ^ oo. The best fit has a slope of —0.9, close to —1 instead of —4/3, which is explained near the end of SectionUl 




FIG. 2: (Color online) Results for the tethered chain. A log-log plot of Same as for Fig. [B but using the continuum heat 
current given by Eq. ([26j . 
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TABLE I: Parameters of the tethered chain simulations for different values of particle number A'' and time step St. tcquii is 
the time for equilibration, fmoas is the subsequent time during which measurements are performed, and rirun, is the number of 
"runs", where one run comprises the equilibration plus the measurement time steps. 
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FIG. 3: (Color online) Results for the tethered chain. A log-log plot of Results for constant pressure. A log-log plot of the 
ratio of the continuum thermal conductivity with the heat current defined in Eq. (|13|) to the lattice thermal conductivity with 
the heat current defined in Eq. p2|l minus one. The best fit has a slope of about —0.8. 
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FIG. 4: (Color online) Results for the tethered chain. A log- log plot of Same as Fig. |3]but using the continuum heat current 
given in Eq. (|26p . The error bars are sufficiently large that it is not possible to estimate the size dependence reliably. 



A. Tethered Chain 

Figure [T] presents data for the tethered chain for the ratio of the thermal conductivity from the continuum current 
in Eq. p3|) to that from the lattice current in Eq. ([12]) minus 1, as a function of system size N, for different values of 
the time step St. If the continuum and lattice thermal conductivities agree in the thermodynamic limit, the results 
should tend to zero for N oo. It is clearly essential to extrapolate the results to small values of 6t, and when one 
does so, the result is close to the expected 1/A^ dependence. The best fit of the data for the smallest value of 5t gives 
with r/ ~ 0.9, but the small difference in the value of the exponent r] from one is probably due to corrections to 
scaling which are not completely negligible for this range of sizes. 

In Fig. [2] we show data similar to Fig. [1] but using the second definition of the continuum current, in Eq. (|26|) . 
The difference in thermal conductivities is smaller, so it is harder to get good statistics on it, but it also seems to 
decay with size like 1/7V. This conclusion disagrees with the naive expectation of Eq. ((22|). but is consistent with the 
physical argument given in the introduction. 

B. Constant Pressure Chain 

For the constant pressure chain. Fig. [3] plots the ratio of the thermal conductivity using the continuum heat current 
in Eq. to that using the lattice current in Eq. minus 1. For reasons that are not clear to us, the dependence 
on time step is much smaller than for the corresponding tethered chain results in Fig. [TJ A fit to the data gives a 
size dependence of N'"^ with rj ~ 0.8. Taken naively, this is a faster decay than the l/N^^^ dependence expected 
for constant pressure, as discussed in the introduction. However, the error bars on the exponent are sufficiently large 
that any 1 > > i is possible. It is also possible that larger sizes are needed to see the asymptotic size dependence. 
We note that the relative uncertainty in the length of the system is of order l/N^^^ and it seems surprising to us that 
difference in thermal conductivities could be smaller than this asymptotically. Despite these uncertainties, the relative 
difference in thermal conductivities is small, even for constant pressure, and appears to vanish in the thermodynamic 
limit. 

Finally, Fig. [4] is similar to Fig. [3] but uses the alternative definition of the head current in Eq. ([26]) . There is still 
a non-zero difference compared with the lattice conductivity but it is smaller, so the error bars are larger, as a result 
of which it is not possible to reliably estimate any functional form for the difference. 
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V. CONCLUSION 

In this paper, we have shown that the continuum and lattice versions of the heat current yield different results for 
the heat conductivity of a finite oscillator chain with Langevin baths through the Green-Kubo formula. The results 
using the continuum and lattice currents are different for a finite chain. Since the thermal conductivity obtained using 
the lattice current is exact @ , this implies that results using the continuum current are approximate even though this 
current is apparently exact. Numerically, if the error is fit to a form ~ 1/A^'' for a chain of N particles, we obtain 
ry < 1 for a chain that is tethered just beyond its end points. We can argue that the error is at least partly because 
the length of the chain fluctuates, and therefore must be bounded below by ~ 1/iV; our numerical results are close 
to and slightly above this lower bound. The same argument would yield a lower bound to the error of ~ 1/ ^/N for 
a chain at constant pressure. Surprisingly, though, the numerical results for such a chain show an error that decays 
faster, almost as 1/N, although the asymptotic large N exponent may approach 0.5. A small change in the definition 
of the continuum heat current, which is benign in the thermodynamic limit, reduces the error substantially. 
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